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Abstract 


In the current study, changes in the amount of extreme precipitation with different return periods for the 
period 2026—2050 were analyzed using the weighted combination of 10 CMIP6 models under different 
scenarios. The spatial and temporal ranking of the historical simulation performance of the 10 CMIP6 
models in Iran was used to determine the weights for each model. Subsequently, a rank-based weighting 
method was employed to project extreme precipitation levels in the study area. 

For extreme precipitation analysis, the best statistical distribution was selected for each station based on 
historical data. In 84% of the stations, the Log-Pearson Type 3 distribution was identified as the best fit 
using the chi-square test for both the historical period (1990—2014) and future scenarios. 

A comparison between the results of extreme precipitation for the historical and future periods revealed 
that the percentage of change in projected extreme precipitation was notably higher under the SSP2-4.5 
scenario compared to the other two scenarios. Moreover, with longer return periods, the percentage increase 
in extreme precipitation was greater, and the spatial extent of its occurrence also expanded. The greatest 
increases in projected extreme precipitation, exceeding 25%, were observed in the SSP1-2.6 scenario in 
Gorgan, Babolsar, Bojnord, Arak, Isfahan, and Zahedan stations. In the SSP2-4.5 scenario, the largest 
increases were recorded in Gorgan, Sanandaj, Hamedan, Arak, Kashan, Isfahan, Shiraz, Fasa, Minab, 
Bushehr, Bandar Abbas, Bandar Lengeh, and Abu Musi stations. Under the SSP5-8.5 scenario, the greatest 
increases occurred in Bojnord, Birjand, Arak, Kashan, and Semnan stations. 

A significant increase in extreme precipitation was consistently observed in Arak station across all three 
scenarios. For the 200-year return period, extreme precipitation at Arak station increased by approximately 


70% under the SSP1-2.6 and SSP2-4.5 scenarios, and by about 50% under the SSP5-8.5 scenario. 


Keywords: Extreme Precipitation Events, Climate Projection Models, Return Period Analysis, CMIP6 
Models, Shared Socioeconomic Pathways (SSPs), Statistical Downscaling. 
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در پژوهش حاضر تغییرات مقدار بارش فرین با دوره بازگشت‌های مختلف در دوره ۲۱۲۶-۲۰۵۰ با ترکیب وزنی ۱۰ 
مدل CMIP6‏ در سناریوهای مختلف مطالعه شد. Gb‏ رتبه‌بندی مکانی و زمانی ظرفیت‌های شبیه‌سازی تاریخی ۱۰ 
مدل CMIP6‏ در ایران» وزن‌های مربوط به هر یک از مدل‌ها به دست آمد. سپس با استفاده از روش وزن‌دهی مبتنی بر 
رتبه» پیش‌نگری مقدار بارش فرین در منطقه موردمطالعه انجام شد. در ادامه بارش‌های فرین در بهترین توزیع آماری 
منطبق بر داده‌ها برای هر یک از ایستگاه‌ها که در ۸۴ درصد ایستگاه‌ها توزیع 3 log-Pearson‏ با آزمون کای- 
اسکوثر انتخاب گردید. در دوره تاریخی (۲۰۱۴-۱۹۹۰) و آینده در هر سه pyle‏ و دوره بازگشت‌های مختلف 
محاسبه شد. مقایسه پیش‌نگری بارش‌های فرین با دوره تاریخی نشان داد که درصد تغییرات مقدار بارش فرین پیش - 
نگری شده در سناریو SSP2-45‏ نسبت به دو سناریو دیگر بیشتر است. همچنین به سمت دوره بازگشت‌های 
طولانی‌تر درصد افزايش مقدار بارش پیش‌نگری شده افزایش و منطقه رخداد OT‏ نیز گسترش می‌یابد. بیشترین افزایش 
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اراک اصفهان و زاهدان در سناریو the‏ در ایستگاه‌های ISS‏ سنندج همدان, اراک» کاشان» اصفهان, شیراز فساء 
میناب. Gedy‏ بندرعباس. بندرلنگه. ابوموسی و در سناریو بدبینانه در ایستگاه‌های بجنورد. بیرجند. اراک کاشان و 
سمنان مشاهده شد. افزایش مقدار بارش فرین بیش از YO‏ درصد در ایستگاه اراک در هر سه سناریو (حدود ۷۰ درصد 
5s‏ ار قاش 98۴1-26 و 89245 dude‏ +۵ ورد 35 cra ilo‏ 9985-85 در درن ار anf‏ اسان 
مشاهده شد. 


کلیدواژه‌ها: رویدادهای بارش شدید.مدل‌های پیش‌بینی اقلیمی.تحلیل دوره باز گشت.مدل‌های 


6 مسیرهای اجتماعی-قتصادی مشترک ریزگردانی آماری. 


۱-مقدمه 


اثر این رویداد طبیعی به خطر می‌افتد. وفوع این رویداد طبیعی» اندازه و تکرار آن ناشی از عوامل متعددی 


است که گرمایش جهانی سبب افزایش رویداد آن در مناطق مختلف جهان شده است. 


درک و پیش‌نگری بارش‌های فرین به‌عنوان یکی از عوامل رخداد سیل برای ارزیابی 388 تغییر اقليم» 
سازگاری و آسیب‌پذیری از اهمیت قابل‌توجهی برخوردار است. مطالعات متعددی گزارش کرده‌اند که 
رویدادهای بارش فرین در قرن گذشته بیشتر شده‌اند و در قرن بیست و یکم نیز بیشتر اتفاق خواهند افتاد 
(برای مثال Peng et al., 2018; Dike, Lin & Ibe, 2020; Lee, Paek, Park & Boo, 2020; IPCC, 2020; Tang,‏ 
Jo (Hu & Duan, 2021‏ اصلی رخداد این بارش‌های فرین افزایش دما است. وقتی دما افزايش eb us‏ 
رطوبت ویژه اشباع P‏ بیشتر می‌شود و بنابراین Gb‏ رابطه کلازیوس -کلاییرون ) Ruckstuhl, Philipona,‏ 
(Morland & Ohmura, 2007‏ هوا می تواند حاوی مقدار بیشتری بخارآب باشد. بر اساس استدلال‌های نظری» 
میزان رطوبت جو به میزان ۷ درصد بر درجه کلوین در گرمایش جهانی بر طبق رابطه کلازیوس-کلاپیرون 
افزايش می‌یابد که منجر به بارش فرین شدیدتر می‌شود ) ;2007 Allen & Ingram, 2002: Trenberth et al,‏ 
(Zhang & Zhou; 2019‏ درواقع هنگامی که شرایط رخداد بارش ایجاد می‌شود. هوای اشباع بیشتر منجر به 


(Berg, Moseley & Haerter, 2013) می‌شود‎ pes بارش‎ 
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مدل‌های اقلیمی جهانی به‌عنوان ابزارهای کارآمدی برای پیش‌نگری تغیبرات آتی در اقلیم متوسط و فرین 
هستند )2013 (Chen, Sun, Chen & Zhou, 2012: Chen,‏ مدل‌های اقلیمی یکی از ابزارهای adsl‏ برای 
محققان برای درک چگونگی تغییر آب‌وهوا در گذشته و در آینده است. این مدل‌ها فیزیک» شیمی و 
زیست‌شناسی cm‏ زمین و اقیانوس‌ها را با جزئیات بسیار شبیه‌سازی می کنند. مدل‌های اقلیمی با توسعه همه 
جانبه‌ای که داشته‌اند اما همچنان عدم قطعیت زیادی دارند و برای کاهش این عدم قطعیت میانگین همادی 
چند مدل برای پیش‌نگری ازجمله پیش‌نگری دما و یا بارش استفاده می‌شود ( Sarabi, Dastorani & Zarrin,‏ 


2020; Sarabi, Dastorani & Zarrin, 2021; Zarrin & Dadashi-Roudbari, 2021a; Zarrin & Dadashi- 

.(Roudbari, 2022a; Zareian, 2022; Khansalari & Mohammadi, 2023‏ 
مدل‌های CMIP6‏ بهترین عملکرد را در تخمین بارش و تغییرات مکانی-زمانی آن در مناطق خشک و بسیار 
خشک کشور Ol ul‏ از خود OUS‏ می‌دهند. باوجوداین ارزیابی عملکرد تخمین‌های بارش از مدل‌های مختلف 
;3CMIPG‏ مناطق اقلیمی مختلف ایران OUS‏ می‌دهد که هر یک از مدل‌های 01۷11۳6 به نظر می‌رسد که برای 
یک منطقه اقلیمی خاص. مناسب‌تر باشد )2021 (Yazdandoost, Moradian, Izadi & Aghakouchak,‏ در 
رابطه با کاربست مدل CMIP6‏ در TIU‏ بارش فرین زرین و داداشی رودباری Zarrin & Dadashi-)‏ 
Roudbari,‏ ۰620222 به بررسی روزهای همراه با بارش em‏ کین در ol ul‏ با استفاده از سه مدل EC-Earth3‏ « 
EC-Earth3-CC‏ و EC-Earth3-Veg‏ با روش میانگین وزنی با رویکرد مستقل ums‏ نتایج این مطالعه 
نشان داد بارش سنگین در ایران تحت شرایط تغییر اقلیم آینده افزايش خواهد داشت. زرین و داداشی 
رودباری (Zarrin & Dadashi-Roudbari, 2021b)‏ به مطالعه شدت بارش شدید در ایران بر اساس مجموعه 
پنج مدل GFDL-ESM4, IPSL-CM6A-LR, MPI-ESMI-2-HR, MRIESM2-0, UKESMI-0-) CMIP6‏ 
(LL‏ پرداخته‌اند و جهت ادغام نتایج حاصل از چندین مدل. از روش وزن‌دهی مستقل استفاده شده است که 
در بین این مدل‌ها مدل MPEESMI-2-HR‏ بالاترین عملکرد را OUS‏ می‌دهد. همچنین زرین و همکاران 
(Zarrin, Dadashi-Roudbari & Hassani, 2022c)‏ بر اساس مجموعه همان پنج مدل CMIP6‏ در coulis‏ پیش = 
نگری می کند که آستانه‌های بارش برای هر دو شاخص ۸10۳۳ و R20mm‏ در دهه‌های آینده در سواحل 


دریای خزر افزایش خواهد یافت Jy‏ شاخحص‌های بارش برای R95p‏ در مناطق غربی و جنوب غربی ; R99p‏ 
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در مناطق جنوبی ایران کاهش خواهند داشت. همچنین تغییرات در شاخص‌های بارش فرین در حوزه دریاچه 
ارومیه با استفاده از ۱۶ مدل CMIPS‏ توسط جوان و همکاران )2023 (Javan, Movaghari & Park,‏ 
موردبررسی قرار گرفت. نتایج این پژوهش نشان دادند که بارش میانگین حوضه تا پایان قرن ۲۱ کاهش و 
شاخص روزهای خشک متوالی افزایش خواهد یافت. بااین حال. سایر شاخص‌های بیشینه بارش یک‌روزه و 
eine‏ روزهای بسیار مرطوب. روزهای متوالی مرطوب و بارش روزهای مرطوب نسبت به دوره مرجع 
کاهش خواهد Cal‏ همچنین خوانسالاری و محمدی )2023 (Khansalari & Mohammadi,‏ پیش‌نگری 
بارش فرین در کشور ايران را با استفاده از پنج مدل EC-Earth3- EC-Earth3 MPI-ESM1-2-HR) CMIP6‏ 
GFDL-ESM4 .Veg‏ و (MRI-ESM2-0‏ موردبررسی قرار داده‌اند. نتایج این مطالعه نشان داد افزایش مقدار 
بارش فرین بر افزايش تعداد روزهای وقوع بارش فرین تقدم دارد؛ به‌بیان‌دیگر این بارش‌های فرین در تعداد 


روزهای کمتری رخ خواهند داد که خود اعلام‌خطری برای رخداد بارش‌های سیل آسا می‌باشد. 


برای پیش‌نگری دوره بازگشت بارش فرین روزانه در ایران از داده‌های مشاهداتی و کاربست یک مدل اقلیمی 
و یا میانگین همادی چند مدل و یا بررسی کارایی پایگاه داده بارش همادی-وزنی چند منبعی (MSWEP:‏ 
Multi-Source Weighted-Ensemble Precipitation)‏ و کاربست آن در تعیین دوره‌های بازگشت POPE‏ 
احتمال رخداد بارش سالانه در ایران نیز مطالعاتی انجام شده است که این مطالعات در مقایسه با بررسی 
بارش فرین بسیار e$‏ است ( ,20220 Hong, Javan, Shin & Park, 2021; Zarrin & Dadashi-Roudbari,‏ 
(Azizi & Nejatian, 2‏ هونگ و همکاران )2021 (Hong et al.,‏ برای پیش‌نگری بارش فرین در ایران از 
داده‌های مشاهداتی و مدل همادی CMIP6‏ استفاده نموده‌اند. نتایج این پژوهش OU‏ می‌دهد که در ایران با 
دوره بازگشت ۲۰ ساله سطح بارش S‏ روزانه در سه سناریوی SSP245‏ 58۳370 و 59۳585 تا سال 
۰ به ترتیب ۱۵/۶ ۲۳/۲ و ۲۸/۷ درصد افزایش می‌پابد. همچنین زرین و داداشی رودباری ) & Zarrin‏ 
(2022b Dadashi-Roudbari,‏ پژوهشی با هدف بررسی کارایی پایگاه داده بارش همادی-وزنی چند منبعی 
(MSWEP: Multi-Source Weighted-Ensemble Precipitation)‏ و Cus jS‏ آن در One‏ دوره‌های بازگشت 


و بررسی احتمال رخداد بارش سالانه در ایران انجام داده‌اند. نتایج این مطالعه نشان می‌دهد ud 4S‏ پایگاه داده 
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در مناطق خشک و نیمه‌حشک کشور بهترین کارایی و در پهنه پرباران و کوهستانی کمترین کارایی را دارد. 
همچنین بیشینه مقدار بارش محتمل در دوره‌های Y‏ تا ۱۵ ساله در سواحل شمالی و auis‏ شدت آن در 
سواحل جنوبی ایران اتفاق می‌افتد. در این مطالعه احتمال رخداد بارش در ایران بین ۲/۰۳ تا ۳۲/۰۲ درصد 
بیان شده است. همچنین عزیزی و نجاتیان )2022 (Azizi & Nejatian,‏ به بررسی شدت- مدت و دوره 
5b‏ گشت شاخحص‌های خشکسالی SPI‏ و SPEI‏ با مدل EC-EARTH‏ در دشت ورامین پرداخته‌اند. نتایج این 
پژوهش نشان می‌دهد که شدت خشکسالی در آینده نسبت به دوره ah‏ بر اساس شاخحص‌های SPI‏ و SPEI‏ 
به ترتیب ۸ و YA‏ درصد افزايش دارد. همچنین دوره بازگشت خشکسالی‌های شدید تحت سناریوی 10۳8.5 
بر اساس شاخحص SPEI‏ در دوره‌های al‏ و آتی به ترتیب ۸ و ۶ سال است که نشان دهنده کاهش دوره 


بازگشت خشکسالی‌های شدید و افزایش سال‌های خشک در آینده است. 


با توجه به تغییر اقلیم و رویه افزایش دمای جهانی» درحالی که بسیاری از مناطق به‌احتمال‌زیاد خشک‌تر می- 
شوند. مطالعه دوره بازگشت و احتمال رخداد بارش‌های فرین اهمیت زیادی دارند؛ اگرچه دوره‌های خشک 
طولانی‌تر این رویدادهای منفرد را از هم جدا می‌کنند )2015 (Mann & Kump,‏ با توجه به افزايش دوره 
خشکی و کوتاه شدن دوره وقوع بارش و رخداد این بارش‌ها به‌صورت فرین نیاز به برنامه‌ریزی و مدیریت 
جهت حفاظت و بهره‌برداری بهینه از آب در دسترس در زمان رخداد فرین‌ها و کاهش خطرات بارش فرین 
dole‏ سیل‌های گسترده است. در این مطالعه با هدف بررسی تغییرات مقداز باوش فرین با دوره باز گشت‌های 
مختلف (۲ ۸۵ ۸۱۰ ۲۵ ۸۵۰ ۱۰۰ و ۲۰۰ سال) در ایران بر اساس رویکرد همادی مدل‌های CMIP6‏ در آینده 
نزدیک» از روش وزن‌دهی نوین برای اختصاص به مدل‌های مختلف به‌کار رفته در پیش‌نگری بارش فرین و 
از شاخص‌های مناسب به‌منظور تعیین بارش‌های فرین در bE‏ مختلف ایران استفاده شده است؛ بنابراین 
هدف اصلی پژوهش حاضر تعیین مقدار بارش فرین در دوره بازگشت‌های مختلف و تغییرات آن نسبت به 
agus‏ تازبخی تا dau‏ قرن ۲۷ با ترکیب وز V‏ مدل مورداستفاده از مدل‌ها ۷1۳6و God‏ ماز رهاق 


SSP2-4.5 SSP1-2.6‏ و SSP5-8.5‏ است. 


۱۳۰۳ جلد ۰۱۳ شماره ۲ پائیز‎ à deoa نشریه جغرافیا و مخاطرات‎ yy. 


۲-مواد و روش‌ها 
۱-۲ منطقه موردمطالعه 


ایران» واقع در منطقه‌ای با اقلیم گذرای بین مناطق خشک و نیمه‌حشک تغییرات قابل‌توجهی در الگوهای 
بارش به‌ویژه از نظر مکانی و زمانی OUS‏ می‌دهد. تنوع کوهساری کشور که شامل bee Sars,‏ فلات‌ها و 
دشت‌ها می‌شود. بر اقلیم محلی و سامانه‌های جوی تأثیر می‌گذارد و منجر به تفاوت‌های بارز در توزیع بارش 
می گردد. مطالعات اخیر بر شدت و فراوانی فزاینده‌ی رویدادهای بارش شدید ASE‏ کرده‌اند که تحت تأثیر 
تغییرات اقلیمی و نوسانات جوی قرار دارند. با توجه به اینکه ایران به‌شدت در برابر تأثیرات شرایط جوی 
شدید. ازجمله سیلاب‌ها و خشکسالی‌ها آسیب‌پذیر است. درک تغییرات مقدار بارش‌های شدید برای 
مدیریت مژثر منابع آب و کاهش خطرات ناشی از بلایای طبیعی ضروری است. این مطالعه با استفاده از یک 
رویکرد همادی و به کارگیری ده مدل۷11۳6) به پیش‌بینی تغییرات در polis‏ بارش شدید با دوره‌های 
بازگشت مختلف می‌پردازد. نتایج این تحقیق به‌منظور ارائه بینش‌هایی در مورد چالش‌های آب‌شناسی آینده 
ایران و کمک به درک بهتر تأثیرات تغییرات اقلیمی بر آب‌شناسی منطقه‌ای و اطلاع‌رسانی به سیاست‌گذاران 


۲-۲ داده‌های مورداستفاده 


با توجه به عدم قطعیت مدل‌ها )2011 (Hawkins & Sutton,‏ برای افزایش سطح قابلیت اطمینان در این 
بررسی از روش همادی چند مدلی (MME!)‏ برای پیش‌نگری اقلیمی استفاده شده است؛ بنابراین در این 
مطالعه از مجموعه داده‌های تاریخی (۲۰۱۴-۱۹۹۰) و پیش‌نگری (۲۰۵۰-۲۰۲۶) بارش ۱۰ مدل از مدل‌های 
۳ استفاده شده است که مشخصات مدل‌ها در جدول Y‏ ارائه شده است. با توسعه همه‌جانبه مدل‌های 
اقلیمی. همچنان این مدل‌ها عدم قطعیت زیادی دارند و برای کاهش این عدم قطعیت میانگین همادی چند 
مدل برای پیش‌نگری متغیرهای هواشناسی استفاده می‌شود؛ بنابراین انتخاب مناسب این مدل‌ها برای منطقه 


1 Multi Model Ensemble 
2 https://esgf-node.IInl.gov/search/cmip6/ 
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موردمطالعه با توجه معیارهای مختلف و افزایش تعداد مدل‌های دخالت داده شده در تحقیق از اهمیت 
بسیاری برخوردار است. در انتخاب مدل‌ها در این مطالعه. به ۳ مورد تفکیک افقی "NV.‏ حساسیت اقلیمی و 


نتایج مطالعات ome‏ توجه شده است. 


۱. بر اساس تعریف تفکیک‌های افقی متنوع در مطالعه )2022 (Lin & Yu,‏ که شامل تفکیک افقی بالا 


Y. 5l yu)‏ کیلومترمربع؛ ۲ مدل)» سعی شده است که تور T‏ از تمام تفکیک‌های افقی استفاده شود. 


۲ حساسیت اقلیمی "Climate Model Sensitivity"‏ معنای تغییر دما در پاسخ به تغییر در ترکیب 
گازهای گلخانه‌ای در جو زمین است. این مفهوم در علوم اقلیمی به‌عنوان "حساسیت دما به افزايش دی‌اکسید 


کربن نیز شناخته می شود. حساسیت اقلیمی معمولاً به دو نحو اندازه‌گیری می‌شود: 


الف: (TCR - Transient Climate Response)‏ این معیار نشان‌دهنده تغییر دما در دوران کوتاه‌مدت است که 


این دوره معمولا به مدت چند دهه تا یک قرن است. 


ب: (ECS - Equilibrium Climate Sensitivity)‏ این شاخحص نشان‌دهنده تغییر دمای مدل‌های اقلیمی به عدم 
تعادل بلندمدت دما در پاسخ به افزایش دی‌اکسید کربن است. درواقع. ECS‏ تنشان‌دهنده تغییر دمای نهایی و 


تابتی است که بعد از دوره بلندمدت اتفاق می‌افتد. 


TCR polis‏ و ECS‏ برای هر مدل اقلیمی میزان حساسیت آن مدل به افزایش دی‌اکسید کربن را نمایش 
می‌دهند. افزایش دی‌اکسید کربن از انشعاب‌های انسانی مانند سوخت‌های فسیلی و فعالیت‌های صنعتی Line‏ 
می گیرد و باعث افزایش دما و تغییرات اقلیمی جهانی می‌شود. حساسیت اقلیمی برای پیش‌بینی تغییرات 
اقلیمی آینده و ارزیابی تأثیر افزایش دی‌اکسید کربن اهمیت دارد. به‌عنوان‌مثال» اگر مدل اقلیمی OUS‏ دهد که 
افزایش دی‌اکسید کربن در جو زمین تغییرات دمایی بسیار بزرگی ایجاد می‌کند. حساسیت اقلیمی این مدل به 
افزایش دی‌اکسید کربن بالا است؛ به‌عبارت‌دیگر مدل اقلیمی در این حالت بسیار حساس به تغییر دی‌اکسید 
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اقلیمی در آینده است. مدل‌هایی که به‌درستی حساسیت‌های اقلیمی را نمایان می کنند می‌توانند پیش‌بینی‌های 


دقیق‌تری ارائه دهند. 


میل و همکاران )2020 (Meehl et al.,‏ به بررسی در مورد حساسیت اقلیم ترازمند (ECS)‏ و پاسخ اقلیمی گذرا 
(TCR)‏ در مدل‌های سیستم زمین CMIP6‏ پرداخته‌اند. مدل‌های مورد ذکر در جدول ۱ در متن مقاله مورد 
داوری» با تفکیک افقی و حساسیت اقلیم ترازمند مختلف (ECS)‏ انتخاب شده‌اند. بر اساس مطالعه پیشین از 
نظر حساسیت Cas (ceil‏ مدل INM-CM5-0 MRI-ESM2-0 GFDL-ESM4 MPI-ESM1-2-HR‏ 
MPI-ESM1.2-LR 6‏ و NorESM2-MM‏ با حساسیت اقلیمی کم (با حساسیت اقلیمی کمتر از 
میانگین ۳۹ مدل موردمطالعه) و سه مدل EC-Earth3-Veg EC-Earth3‏ و CanESMS‏ با حساسیت اقلیمی 
365 (با حساسیت اقلیمی بیشتر از میانگین YA‏ مدل موردمطالعه) انتخاب مناسبی برای کل مجموعه CMIP6‏ 
هستند. با توجه به حجم زیاد محاسبات سعی شده است انتخاب مناسبی از مدل‌ها با تفکیک افقی و 


حساسیت اقلیمی مختلف انجام شود. 


Tang et al., 2021; Zarrin & Dadashi-) „>| در انتخاب مدل‌های مذکور همچنین به نتایج مطالعات‎ Y 


Roudbari, 2021b; Zarrin & Dadashi-Roudbari, 2022a; Zarrin et al, 2022c; Khansalari & 
نیز توجه شده است.‎ (Mohammadi, 2023 


همچنین به‌منظور ارزیابی و تعیین وزن هر یک از این مدل‌ها در این مطالعه از مجموعه داده‌های بارش 
دیدبانی شده و داده‌های تاریخی مدل‌های مورداستفاده CMIP6‏ در بازه زمانی مشترک (۲۰۱۳-۱۹۹۰) استفاده 
شده است. منطقه موردمطالعه (شکل ۱) کشور ایران ably‏ در جنوب‌غرب آسیا و تقریباً بین ۲۶ تا ۴۰ درجه 
شمالی و ۴۴ تا ۶۲ درجه طول جغرافیایی می‌باشد. موقعیت ایستگاه‌های همدیدی موردمطالعه برای تعیین 
دوره بازگشت و احتمال رخداد بارش‌های فرین تا اواسط قرن ۲۱ با ترکیب وزنی ۱۰ مدل مورداستفاده از 
مدل‌های CMIP6‏ و تحت سناریوهای SSP2-4.5 SSP1-2.6‏ و SSPS-8.5‏ و کوهساری منطقه در شکل ۱ 


OUS‏ داده شده است. انتخاب ایستگاه‌ها بر اساس وجود دوره آماری داده‌های دیدبانی ایستگاه‌ها در دوره 


yyy 
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تاریخی مدل‌های CMIP6‏ (۱۹۹۰-۲۰۱۴) و حداکثر pue‏ وجود داده در ایستگاه مورد نظر برابر با ۰ درصد 


40N 


35N 


30N 


25N 


صورت گرفته است. 


Height (m) 


50E 55E 60E 


شکل ۱- موقعیت ایستگاه‌های همدیدی موردمطالعه به همراه کوهساری منطقه 


Fig.1. Location of the synoptic stations studied, along with the mountainous areas of the 


region 


حدول ۱- مشخصات مدل‌های CMIP6‏ مورداستفاده 


Table 1 - Specifications of the CMIP6 models used 


مو سسه سازنده 


) (تفکیک افقی و تعداد تراز قائم‎ resolution ردیف نام مدل‎ 
(Manufacturing company) and number of vertical levels) (Model name) (Row) 
. 0.9375 #0.9375 
مرکز داده‌های اقلیمی المان‎ 
تراز: ۰/۰۱ هکتویاسکال‎ op SVL Gls ۵ 
German Climate Data Center iid € xod : MPI-ESM1-2-HR 1 
(Deutsches Klimarechenzentrum) 95 Elevation, The highest elevation: 
0.01Hpa 
اتحادیه‌ای از چندین شرکت تجاری‎ 0.7031 #0.7031 
کشور اروپایی)‎ Y) تراز بالاترین تراز: ۰/۰۱ هکتوپاسکال‎ ۱ EC-Earth3 " 
Consortium of Multiple s . ۲ 
: ۱ 91 Elevation, The highest elevation: 
Commercial Companies (12 0.01H 
European Countries) ini ia 
اتحادیه‌ای از چندین شرکت تجاری (۱۲ کشور‎ 0.7031 «0.7031 
تراز بالاترین تراز: ۰/۰۱ هکتویاسکال اروپایی)‎ ۱ EC-Earth3-Veg 3 


Consortium of Multiple 
Commercial Companies (12 


9] Elevation, The highest elevation: 
0.01Hpa 
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European Countries) 


اداره ملی اقیانوسی و جوی» آزمایشگاه 
دینامیک سیالات ژئوفیزیکی آمریکا 
National Oceanic and Atmospheric‏ 
Administration (NOAA),‏ 


Geophysical Fluid Dynamics 
Laboratory (GFDL), USA 


موسسه تحقیقات هواشناسی ژاپن 
Japan Meteorological Research‏ 
Institute (MRI)‏ 


مرکز کانادایی مدل‌سازی و تحلیا اقلیم. 
Canadian Centre for Climate‏ 
Modelling and Analysis (CCCma),‏ 
Environment and Climate Change‏ 
Canada‏ 
موسسه ریاضیات عددی, آکادمی علوم روسیه 
Institute of Numerical‏ 
Mathematics, Russian Academy of‏ 
Sciences (INM RAS)‏ 

علوم و فناوری‌های دریایی-زمینی آژانس «qu‏ 
موسسه تحفیقات جوی و اقیانوسی دانشگاه 
توکیو ژاپن» موسسه ملی مطالعات 
زیست‌محیطی ژاپن و مرکز RIKEN‏ برای 
علوم محاسباتی ژاپن 

Japan Agency for Marine-Earth 
Science and Technology 
(JAMSTEC), Atmosphere and 
Ocean Research Institute at the 
University of Tokyo, National 
Institute for Environmental Studies 
(NIES), and RIKEN Center for 
Computational Science, Japan 


موسسه هواشناسی ماکس پلانک آلمان 


Max Planck Institute for 
Meteorology of Germany 


اتحادیه مدل‌سازی اقلیم نروژ 


Norwegian Climate Modelling 
Consortium 


IE 
تران بالاترین تراز: ۰/۰۱ هکتوپاسکال‎ Y4 


49 Elevation, The highest elevation: 
0.01Hpa 


1.125% 1.125 
بالاترین تراز: ۰/۰۱ هکتوپاسکال‎ GL ۰ 
80 Elevation, The highest elevation: 
0.01Hpa 
2.8125 5 
هکتوپاسکال‎ Y تران بالاترین تراز:‎ ¥4 


49 Elevation, The highest elevation: 
1Hpa 


2« 1.5 
۳ تران بالاترین تراز سیگما: ۰/۰۰۰۲ 


73 Elevation, The highest Sigma: 0.0002 


1.4063« 1.4063 
GL ۱‏ بالاترین تراز: ۰/۰۰۴ کتوپاسکال 


81 Elevation, The highest elevation: 
0.004 Kpa 


1.875% 1.875 
بالاترین تراز: ۰/۰۱ هکتوپاسکال‎ GL ۷ 
47 Elevation, The highest elevation: 
0.01Hpa 
0.9375 « 1.25 
هکتوپاسکال‎ Y تراز بالاترین تراز:‎ ۲ 
32 Elevation, The highest elevation: 3 
Hpa 


GFDL-ESM4 


MRI-ESM2-0 


CanESM5 


INM-CM5-0 


MIROC6 8 


MPI-ESM1.2-LR 9 


NorESM2-MM 10 
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۳-۲ روش مطالعه 
الف: محاسبه شاخص‌های بارش فرین 


میانگین بارش ایران به دلیل موقعیت جغرافیایی بسیار کمتر از میانگین جهانی است. میزان بارش ۳۰۰ میلی 
متر در مناطق Kis‏ به‌عنوان حداقل بارش برای کشت ديم عنوان شده است. بیش از دو سوم مساحت ایران 
دارای میانگین بارش سالانه کمثر از ۳۰۰ میلی‌مثر است» اما هسته‌های بارش بالای ۱۰۰۰ میلی‌مثر در منطقه 
زاگرس و دامنه‌های شمالی البرز نیز وجود دارد )2008 (Alijani, O'Brien & Yarnal,‏ با توجه به دامنه وسیع 
میانگین بارش در کشور انتخاب یک آستانه بارش فرین در کل کشور سوال‌برانگیز است؛ بنابراین در این 
مطالعه آستانه بارش فرین با استفاده از صدک‌ها به روش )2005 (Zhai, Zhang, Wan & Pan,‏ محاسبه می - 
شود. بر اساس این آستانه بارش فرین» چهار شاخحص بارش فرین (جدول (Y‏ که توسط کارگروه ETCCDI!‏ 
ارائه و توسط سازمان جهانی هواشناسی کمیسیون اقلیم‌شناسی و توصیه‌های برنامه تغییر اقلیم و پیش‌بینی 
توسعه adh‏ است» مورداستفاده قرار می‌گیرد. این چهار شاخص فرین AEPI 1950 RISP‏ و 9501 به 
ترتیب مجموع بارش فرین بر حسب واحد میلی‌متر تعداد روزهایی که در سال موردبررسی مقدار بارش در 
bol‏ از آستانه بارش فرین Of‏ سال بیشتر باشد. شدت مطلق بارش فرین بر حسب میلی‌متر بر روز (نسبت 
شاخص اول به شاخص دوم) و نسبت بارش ce‏ به بارش سالانه در روزهای بارانی (بارش روزانه بیش از 


یک میلی‌متر) بر حسب در صد. متغیرهای آماری اساسی هستند که بارش فرین را توصیف می کنند. 


حدول ۲- معرفی مختصر جهار شاخص بارش فرین در این مطالعه )2021 (Tang et al.,‏ 


Table 2 - Brief introduction of the four extreme precipitation indices used in this study (Tang et al., 


2021). 
(Index) شاخص‎ (Descriptive Name of Index) نام توصیفی شاخص‎ (Unit) واحد‎ 
R95p Precipitation due to very wet days mm 
R95d Extreme precipitation days day 


! Expert Team on Climate Change Detection and Indices 
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AEPI Absolute intensity extreme precipitation Index mm/d 


R95pT Extreme precipitation to total precipitation of very wet days / 


| فضا در‎ Sleds ah il cl OMS, کد‎ etl dta Ve e Stasi als 
و نسبت انحراف معیار‎ (RMS) از طرح تیلور با بررسی همبستگی الگو» تفاوت ریشه میانگین مربعات‎ 
از اماز‎ Saye خیرات زماتی‎ duos ely 6۳6 cadia کی‎ al بو برای‎ (Taylor 2001) 
مطابق رابطه‎ (Chen, Jiang & Li, 2011) مهارت تغییریذیری بین‌سالانه توصیف‌شده توسط چن و همکاران‎ 


)۱( استفاده می‌شود. 
رابطه (۱) 


STDm STD, 
IVS = (——- 
STD, STD; 


?) 
که در این رابطه STDm‏ و STDo‏ به ترتیب انحراف معیارهای بین‌سالانه شبیه‌سازی‌های مدل و مشاهدات 
هستند. 


ب: رتبه‌بندی و وزن‌دهی هر مدل مورداستفاده از CMIP6‏ 


با توجه به قابلیت‌های شبیه‌سازی مکانی (نمودار تیلور) و زمانی (مقادیر (IVS‏ هر یک از شاخص‌های 
بارش فرین برای رتبه‌بندی مدل‌های CMIP6‏ استفاده می‌شوند ( & Tang et al, 2021: Khansalari‏ 
Mohammadi, 2023‏ بر این اساس» سنجه رتبه‌بندی کل Giang, Li, Xu & Li, 2015) MR‏ به‌صورت رابطه 


(۲) محاسبه می‌شود: 


11 
1 
MR=1 -=> ranki (Y) 
nm بط‎ 
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کر Gil‏ سول هراد Bie eet‏ ار کرت A‏ لاه a Velde‏ عم تنعل براق 1 امین فا 
S ee‏ اسر Gib Rie ee o6‏ مر EG‏ وین Gane)‏ مان snc pb Sy oly n‏ ا 
مقدار MR‏ به ۱ نزدیک‌تر باشد, مهارت شبیه‌سازی بالاتر است. برای هر مدل 0041۳6 دو مقدار MR‏ (یکی با 
توجه به نمودار تیلور و دیگری بر اساس مقادیر IVS‏ به دست می‌آید. سپس برای هر مدل. Si‏ محاسبه می - 
شود که مجموع رتبه‌های دو مقدار MR‏ است. ازاین‌رو وزن متناظر هر مدل به روش چن و همکاران 


Chen)‏ 06 15 201( در رابطه (Y)‏ محاسبه می‌شود: 


M 
R, = 2S ۳ 
زگ‎ 
R; 
= (۴) 
ini Ri 


را دارا هستند. با استفاده CRi 3b‏ می‌توان وزن هر مدل (Wi)‏ را با استفاده از معادله ۴ به دست آورد که درواقع 


W;‏ نرمال شده Ri‏ است. 


TI UNUS تغییرات پیش‌بینی شده» برای اندازه گیری عدم قطعیت‎ (RMSD') cà میانگین مربعات‎ Ane) 
وزن متناظر با آن (محاسبه شده با‎ W مقدار پیش‌بینی مدل ام و‎ P اگر‎ (Chen et al., 2011) استفاده می‌شود‎ 


استفاده از فرمول ۴ باشد. عدم قطعیت پیش‌نگری از روابط ۵ و ۶محاسبه می‌شود: 


B (0) 
iW; 


5 نی‎ Wii - WME)? (۶) 
RMSD = TTIW 


! Root-Mean-Square Difference 


WMME — 


۱۳۰۳ جلد ۰۱۳ شماره ۲ پائیز‎ à deoa نشریه جغرافیا و مخاطرات‎ YYA 


Td‏ دوره بازگشت 


حداکثر بارش ۲۴ ساعته یکی از عناصر اساسی در مطالعات سیلاب‌ها و همچنین طراحی تاسیسات آبی می = 
باشد. بر این اساس» برای محاسبه آن در دوره بازگشت‌های مختلف انتخاب توزیع آماری مناسب ضروری 
است. در بررسی‌های اقلیمی سعی می‌شود داده‌های اندازه‌گیری و ثبت شده با توابع توزیم تثوری برازش داده 
شود تا بهترین تابعی که با داده‌ها مطابقت دارد به‌عنوان تابع توزیع احتمال برگزیده شود )2015 (Alizadeh,‏ 
در این پژوهش برای دوره آماری ۵ ساله (دوره تاریخی ۱۹۹۰-۲۰۱۴ و دوره آتی ۲۰ ۲۰۲۶) در 
ایستگاه‌های دیدبانی کشور ایران» مجموعه آمار به‌دست‌آمده از ترکیب وزنی مبتنی بر رتبه ۱۰ مدل از 
6 به کمک نرم‌افزار تحلیل فراوانی آب‌شناسی Easy Fit‏ با چندین نوع توزیع آماری شامل توزیع‌های 
Gamma Lognormal Lognormal(3p) Normal Weibull Exponential‏ و Log-Pearson type3‏ به‌منظور 
برآورد بهترین توزیع» برازش داده شد. همچنین ارزیابی معنی‌داری توزیع‌های به دست آمده با استفاده از 
آزمون کای اسکوثر (Chi-Squared)‏ انجام شد ) Alijani & Afshar Manesh, 2015; Asadi & Akbari‏ 
(Azirani, 2021‏ در مطالعات مختلف با استفاده از آزمون‌های نکویی. میزان انطباق تابع توزیع با داده‌های 
مشاهداتی تعبین می‌گردد. اگر میزان تطابق در حد قابل قبولی بود. آن توزیع به‌عنوان توزیع برگزیده پذیرفته 
می‌شود .لازم به ذکر است که هیچ توزیع آماری نمی‌تواند دقیقاً بر داده‌های مشاهده‌ای برازش خوبی داشته 
باشد و انتخاب یکی به‌عنوان بهترین توزیع با توجه به مقایسه نسبت به ple‏ توزیع‌ها است ) & Shahabfar‏ 


(Qiyami, 01 


رابطه آزمون CIS‏ اسکوثر به‌صورت رابطه (V)‏ می‌باشد: 
T » (Pei = Poi)? (V)‏ 
Dei‏ 


که ۸ تابعی از درجه آزادی و درصد اطمینان و Dei‏ و Poi‏ به ترتیب مقدار بارش تخمین زده شده و 


مشاهداتی می‌باشند. لازم به ذکر است در این پژوهش درصد اطمینان ۹۵ / در نظر گرفته شد. برای محاسبه 
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دوره بازگشت بارش‌هاء از نرم‌افزار Hyfran-plus‏ استفاده شد. دوره بازگشت به برآورد احتمال روی‌دادن یک 
پدیده مانند سیل یا تغییرات دبی رودخانه گفته می‌شود. معمولاً برآورد دوره بازگشت بر ah‏ اندازه‌گیری 
آماری داده‌های تاریخی به‌منظور دست‌یابی به میانگین زمان تکرار پدیده در یک دوره زمانی است. دوره 
با گشت تعداد سال‌هایی است که به‌طور متوسط بین وقوع دو حادثه مشابه وجود دارد. اگر دوره بازگشت T‏ 


(سال) و احتمال وقوع P‏ باشد. رابطه (A)‏ است )2015 (Alizadeh,‏ 
(A)‏ 


e 
(P 
و یا احتمال ۹۸ درصد یعنی‎ OLS بنابراین احتمال ۱۰ درصد یعنی دوره بازگشت ۱۰ سال (هر ۱۰ سال‎ 


دوره بازگشت ۱/۳۲ üla‏ (در هر ۱۰۲ سال ۰ بار). 
۳- نتایج و بحث 
۱-۳- توزیع مکانی میانگین چهار شاخص بارش فرین در دوره تاریخی ۱۹۹۰-۲۰۱۴ 


شکل Y‏ توزیع مکانی میانگین سالانه چهار شاخص بارش فرین در ایران را که در بخش روش کار به آن‌ها 
اشاره شد. طی دوره تاریخی ۰ G‏ ۲۰۱۴ نشان می‌دهد. بر اساس شکل Ya)‏ شاخحص R95p‏ همراه با 
کوهساری ایران تغییر می‌کند. بیشترین مقادیر شاخحص R95p‏ در سواحل دریای خزر به‌ویژه نواحی غربی و 
جنوبی مشاهده می‌شود که بارش فرین بیش از ۴۰۰ میلی‌متر است. همچنین میزان بارش فرین در دامنه‌های 
غربی رشته‌کوه زاگرس و سواحل شرقی دریای خزر بیش از ۱۰۰ میلی‌متر است. در حقیقت بیشترین رخداد 
بارش فرین در بادسوی دو رشته‌کوه مهم ایران (البرز و زاگرس) رخ می‌دهد. در بادسوی این رشته‌کوه‌های 
مهم در اران که محل ورود سامانه‌های جوی به کشور است با صعود بسته هوا در شیب این رشته کوه‌ها 
فرآیند سرد شدن و اشباع بهتر رخ می‌دهد و منجر به تشکیل ابرهای همرفتی و رخداد بارش بیشتر می گردد. 


همچنین طبق نتایج زرین 5 داداشی رودباری L (Zarrin & Dadashi-Roudbari, 2022b)‏ کاهش عرض 


۳۳۰ نشریه جغرافیا و مخاطرات محیطی, جلد ۰۱۳ شماره ۲ پائیز ۱۳۰۳ 


جغرافیایی از شمال به جنوب و کاهش ارتفاعات از غرب به شرق رخداد بارش کاهش می‌یابد. در شکل b)‏ 
ROS quet sis Y‏ 65 دای اش بر این اسای ما که سا قاد balaa guy‏ بارش 
«cp 5‏ در دامنه‌های شمالی رشته‌کوه البرز و شمال غرب و مناطقی از غرب ایران بیشتر از Y‏ روز است. 
بیشترین Son‏ این شاخص (بیش از ۵ روز) در امتداد سواحل غربی و جنوبی دریای خزر در دوره تاربخی 
۲۰۱۴-۰ مشاهده شده است. بیشترین مقدار شدت بارش فرین نیز در این دوره آماری در سواحل جنوبی 
و غربی دریای خزر و پس از آن در منطقی از جنوب و جنوب غرب کشور (بیش از ۵۰ میلی‌متر بر روز) 
حاصل شده است (شکل (Ve‏ همچنین نسبت بارش فرین به بارش سالانه در روزهای بارانی (بارش روزانه 
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شکل -Y‏ توزیع فضایی میانگین سالانه شاخص‌های بارش فرین برای دوره تاریخی: R95p (a ١۱۹4-۴‏ 


R95pT (d .AEPI (c 18950 (b 


Fig. 2. Spatial distribution of the annual average of the extreme precipitation indices for the 
historical period: 1990-2014. a) R95p, b) R95d, c) AEPI, d) R95pT 


Y-Y‏ - وزن‌دهی مدل‌های مورداستفاده CMIP6‏ بر اساس رتبه‌بندی مکانی و زمانی 


در جدول Y‏ رتبه‌بندی مکانی ۱۰ مدل مورداستفاده 0۷11۳6 برای هر یک از چهار شاخص بارش فرین با 
توجه به سه جنبه نمودار تیلور شامل همبستگی gS‏ نسبت انحراف استاندارد فضایی و ریشه میانگین 
مربعات اختلاف نرمال شده متمرکز و همچنین رتبه‌بندی زمانی بر اساس امتیاز TVS‏ برای اندازه‌گیری مقدار 


شباهت تغییرپذیری بین‌سالانه ارائه شده است. 


جدول ۳- مقادیر آماره‌های تیلور و امتیاز DVS‏ برای هر شاخص بارش فرین» رتبه‌ها و وزن‌های مربوط به ۱۰ 
مدل مورداستفاده CMIP6‏ (جدول )( با توجه به عملکرد جامع آن‌ها با در نظر گرفتن نمودار تیلور و 
تغییرپذ بری بین‌سالانه در کشور ایران برای دوره تاریخی ۱۹۹۰-۲۰۱۴ خانه‌های فرمورنگ شان‌دهنده مدل با 


oer orn 


وزن بزرگ‌تر از ارزش مدل‌های هم‌وزن (۰/۱) است. 
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Table 3 - Values of the Taylor statistics and IVS scores for each extreme precipitation index, 
along with the rankings and weights corresponding to the 10 CMIP6 models used (Table 1), 
based on their overall performance considering the Taylor diagram and interannual 
variability in Iran for the historical period 1990-2014. Red cells indicate models with weights 
greater than the value of equally weighted models (0.1). 


10 9 8 7 6 5 4 3 2 1 Model number 


0.4069 0.3848 0.5644 0.4261 0.4010 0.6390 0.6097 0.8388 0.8742 0.5026 R95p STD ratio 
1.4213 1.29086 | 1.8063 1.7281 1.7199 1.6546 1.5942 1.9598 2.0365 1.2541 R95d STD ratio 
0.3693 0.3075 | 0.3508 0.2025 0.2044 0.4134 0.3741 0.4218 0.4543 0.4999 AEPI STD ratio 
0.8873 1.0892 0.7506 0.7247 1.0701 0.7304 4 0.9596 0.9721 1.1578 — R95pT STD ratio 
1.0092 1.0001 1.1220 1.0247 1.0166 1.1270 1.0079 0.9246 0.9188 0.9468 R95p CNRMSE 
1.5181 1.4202 1.9085 1.8187 1.8296 1.7516 4 1.8516 1.9017 1.3355 R95d CNRMSE 
0.9891 1.0039 1.0200 1.0065 1.0366 0.9584 0.9168 0.8128 0.8082 0.9493 AEPI CNRMSE 
1.3425 1.4598 1.2608 1.2484  Á 1.5474 1.2457 1.2752 1.4107 1.4038 1.4868  R95pT CNRMSE 


0.1799 0.1679 0.0522 0.1535 0.1578 0.1075 0.2913 0.5056 0.5259 0.3538 R95p R 
0.2511 0.2573 0.1711 0.1957 0.1768 | 0.2017 0.2836 0.3597 | 0.3753 | 03141 R95d R 
0.2343 0.2189 0.1113 0.0863 0.0616 0.3190 0.4128 0.6299 0.6361 0.3872 AEPI R 
-0.0091 0.0235 0.0182 0.0238 0.1232 0.0132 0.0221 0.0368 0.0158 4 R95pT R 
2.85 3.09 1.76 4.19 2.95 1.32 1.19 0.71 0.87 1.37 R95p IVS 
1.09 0.57 1.86 2.04 1.29 1.42 0.90 1.16 1.30 0.53 R95d IVS 
11.84 6.35 5.10 10.67 13.58 4.11 4.27 7.52 6.48 2.44 AEPI IVS 
0.27 0.35 0.57 0.90 0.36 0.82 0.39 0.39 0.46 0.30 R95pT IVS 
0.6000 0.5750 0.1250 0.1500 0.0750 0.4500 0.7000 0.5770 | 0.7500 0.6750 MR Taylor 
0.4750 | 0.5250 0.3250 0.0500 0.3000 0.4250 0.6500 0.5500 0.4500 0.7750 MR IVS 
4 5 9 8 10 7 2 5 1 3 Rank of Taylor 
5 4 8 10 9 7 2 3 6 1 Rank of IVS 
9 9 17 18 19 14 4 8 7 4 5 
12.1111 12.1111 64118 6.0556 7.7857 27.2500 13.6250 15.5714 0 R 
0.0904 0.0904 0.0479 0.0452 0.0428 0.0581 0.2035 0.1017 0.1163 0.2035 WwW 


در شکل Y‏ نمودار تیلور چهار شاحص فرین برای ۱۰ مدل مورداستفاده نشان داده شده است. در این نمودار 
مقادیر ضریب همبستگی با امتداد شعاع دایره روی قوس ربع دایره مقادیر انحراف معیار بر روی شعاع ربع 
دایره و مقادیر RMSD‏ بر روی شعاع دوایر متحدالمرکز با مرکز دایره JES‏ مشخص می‌شوند. همچنین 
polis‏ دقیق هر آماره در این جدول oala OLE‏ شده است. با توجه به MR. polis‏ فضایی و زمانی هرچه مقدار 
MR‏ 4 ۱ نزدیی‌تر باشد نشان‌دهنده مهارت بیشتر شبیه‌سازی مدل برای شاخحص‌های بارش فرین است. بر 
این اساس مدل EC-Earth3‏ با MR taylor‏ برابر با ۰/۷۵ (جدول ۳) و با تفکیک افقی ۰/۷۰۳۱۰:۰/۷۰۳۱۰ و 
۱ تراز (جدول ۱) بهترین مهارت در شبیه سازی الگوی فضایی شاخص‌های بارش فرین را دارد که مدل 
مذکور بیشترین تفکیک افقی شبیه‌سازی را در بین ۱۰ مدل مورداستفاده داراست. از دیدگاه شبیه‌سازی زمانی 


(* با ۰/۷۷۵۰ و ۰/۶۵ (جدول‎ ply MR IVS ترتیب با‎ 4 GFDL-ESM4 و‎ MPLESMI2-HR J. 
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مناسب‌تر از ple‏ مدل‌های موردمطالعه هستند. زرین و داداشی رودباری ) Zarrin & Dadashi-Roudbari,‏ 


0 نیز در پیش‌نگری بارش فرین با پنج مدل 011۳6 که سه مدل آن مشترک با پژوهش حاضر است 


مدل MPI-ESMI-2-HR‏ را بهترین مدل معرفی نموده‌اند که در تطابق با نتایج مطالعه حاضر می‌باشد. 


MPI-ESM1-2-HR ® 5 


۰ è MPI-ESM1-2-HR * 5 
* EC-Earth3 ۰ INM-CM5-0 » EC-Earth3 © INM-CM5-0 
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شکل ۳- نمودار تیلور چهار شاخص بارش فر ین برای ۰ مدل مورداستفاده 6 برای دوره تار یخی 


(c) 
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Fig. 3. Taylor diagram of the four extreme precipitation indices for the 10 CMIP6 models used 
for the historical period 1990-2014. a) R95p, b) R95d, c) AEPI, d) R95pT 


پس از محاسبه رتبه ۱۰ مدل CMIPG‏ توسط MR‏ ها بر اساس نمودار تیلور و 1۷5 سپس وزن متناظر lao!‏ 
با استفاده از رابطه ۴ محاسبه شد. بر این اساس وزن مربوط به هر مدل در جدول ۳ ارائه شده است که در 
of‏ خانه‌های قرمزرنگ نشان می‌دهد که وزن مدل‌های مشخص شده بزرگ‌تر از مقدار وزن مدل‌ها در حالتی 
است که وزن هر یک از ۱۰ مدل برابر در نظر گرفته می‌شود (۰/۱). به بیان دیگر مدل‌های مشخص شده تأثیر 


بیشتری در پیش‌نگری احتمالی آینده خواهند داشت, 
-Y-Y‏ نتایج محاسبه دوره بازگشت بارش‌های فرین 
با استفاده از آزمون PL‏ برازش کای اسکوئر تابع (QS‏ احتمالاتی بیشینه بارش ۴ ساعته در دوره‌ی 
تاریخی و دوره ST‏ با سناریوهای مختلف تغییر اقلیم (بر اساس کمترین خطای آزمون)؛ توزیع‌های ارائه شده 


در جدول Y‏ برازش بهتری در مقایسه با سایر توزیع‌های مورد آزمون در ایستگاه‌های موردمطالعه را نشان داد 


و به‌عنوان بهترین توزیع آماری منطبق بر داده‌ها در ایستگاه‌های مورد نظر برای ادامه کار انتخاب شد. 
جدول ۴- توزیع‌های منتخب در ایستگاه‌های موردمطالعه 


Table 4 - Selected distributions at the studied stations 


تابع توزیع ایستگاه 


Station Distribution Function 


Log Pearson type 3‏ اردبیل» Ged‏ ارومیه» رشت» زنجان» سنندج» قزوین» همدان» کرمانشاه» بابلسر GS‏ تهران» اراک» خرم‌آباه 
OU S‏ سمنان» اصفهان» کاشان» شهرکرد. د وگنبدان» ol gal‏ مشهد. یزد. کرمان» سیرجان» بم» شیراز. بوشهرء 


gamma‏ بجنورد» بیرجند» فسا 
di Exponential‏ رلنگه pS‏ 
Weibull‏ 


ابوموسی 


با توجه به تعداد زیاد نمودارها و جداول. درجه مناسب بودن و شاخحص آماری آزمون GIS‏ اسکوثر در 


ایستگاه اردبیل Ol peas‏ نمونه از تابع توزیع با پیشترین عضو )3 «(Log-Pearson type‏ در دوره تاریخی NEP‏ 


۳۳۵ 
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۳ در این جدول» توزیع لوگ پیرسون نوع‎ S سناریو) در حدول ۵ ارائه شده است. بر اساس نتایج‎ Aw) 


برای ایستگاه اردبیل در دوره تاریخی (۲۰۱۴-۱۹۹۰) و آینده (۲۰۵۰-۲۰۲۶) برای هر سه سناریو به دلیل 


کمترین مقدار خطا و مجموع کمترین رتبه انتخاب شد. در جدول ۵ ستون اول انواع توابع توزیم را نشان 


می‌دهد که در مطالعات آب‌شناسی به‌صورت‌های گوناگون با برازش داده‌ها بر آن‌ها مورداستفاده قرار می - 


گیرند )2015 Alizadeh,‏ درحه مناسب بودن. رنه هر توزیع را پس از برازش بر روی داده‌ها oU‏ می دهد 3 


شاخحص آماری» میزان خطا در برآورد توزیع را QUO‏ می‌دهد که هرچه کوچک‌تر باشد مقدار خطا کمتر 


خواهد بود. بر این اساس با توجه به جدول ۱ طبق polis‏ شاحص آماری» کمترین خطا در دوره تاریخی و 


در دوره آینده و در هر سه سناریو مربوط به pU‏ توزیع 3 Log-Pearson‏ است. 


حدول ۵- درحه مناسب بودن و شاخص آماری آزمون کای اسکوئر در ایستگاه اردبیل در دوره تاربخی و آینده 


(سه سناریو) 


Table 5 - Goodness of fit and Chi-square statistical index for the Ardabil station during the 
historical and future periods (three scenarios). 


دوره آینده 2026-2050 


Future period 


دوره تاریخی 
1990-2014 


SSP5-8.5 SSP2-4.5 SSP1-2.6 Historical period 
تلع توزیع درجه درجه درجه درجه‎ 
شاخص شاخص شاخص شاخص‎ Distribution 
t gatti 1 ناسپ‎ n aes i اب‎ Function 
اماری اماری اماری‎ Solel 
Statistica Mid Statistica D Statistica i Statistica bbs 
1 Index Goodnes 1 Index Goodnes 1 Index Goodnes 1 Index Goodnes 
s of Fit s of Fit s of Fit s of Fit 
3.67 1 2.52 1 1.71 1 5.03 1 Log-Pearson 3 
4.13 2 3.12 2 2.36 2 7.51 2 Lognormal (3P) 
8.13 4 10.81 4 6.02 4 7.97 3 Lognormal 
7.55 3 8.98 3 3.34 3 11.09 4 Gamma (3P) 
22.55 7 27.1 7 17.23 7 12.03 5 Weibull 
18.60 5 25.48 6 6.92 5 17.11 6 Weibull (3P) 
21.26 6 20.87 5 13.92 6 28.20 7 Gamma 
23.61 8 31.32 8 24.29 8 41.03 8 Normal 
29.72 9 54.83 9 26.17 9 59.03 9 SE BE 
185.72 10 299.94 10 221.43 10 232.78 10 Exponential 
در این مرحله از مطالعه با توجه به تعریف دوره بازگشت که فاصله زمانی بین رویدادهای با بزرگی یکسان‎ 


است؛ بارش‌های فرین در توزیع‌های منتخب با استفاده از نرم‌افزار Hyfran-plus‏ در هر ایستگاه و در دوره‌ی 


۳۳۶ نشریه جغرافیا و مخاطرات محیطی, جلد ۰۱۳ شماره ۲ پائیز ۱۳۰۳ 


تاریخی و دوره آینده در هر سه سناریو SSP2-4.5 SSP1-2.6‏ و SSP5-8.5‏ با دوره jb‏ گشت‌های مختلف v)‏ 
Yrs Or ۲۵ ۸۰ ۵‏ و ۲۰۰ ساله) محاسبه شد. با توجه ay‏ هدف اصلی این پژوهش و همچنین انبوه 
اطلاعات و توزیع‌های مختلف مقدار بارش فرین پیش‌نگری شده نسبت به دوره تاریخی در شکل X‏ درصد 
تغییرات مقدار بارش فرین پیش‌نگری شده نسبت به مقدار بارش فرین در دوره تاریخی در هر یک از 
ایستگاه‌های موردمطالعه در هر سه سناریو و همه دوره بازگشت‌ها ترسیم شده است. همچنین در این «JS‏ 
در هر دوره و سناریو که مقدار این کمیت بیشتر از YO‏ درصد باشد. ایستگاه مورد نظر با خط قرمزرنگ 
مشخص شده است و در شکل ۵ نمودار تغییرات این کمیت در دوره بازگشت‌های مختلف برای ایستگاه 
مشخص شده ارائه شده است. با توجه به شکل Y‏ در سناریو خوش‌بینانه در سواحل شرقی دریای خزر و 
ایستگاه بجنورد در شمال‌شرق کشور ایستگاه‌های اراک و اصفهان از مرکز کشور و همچنین ایستگاه زاهدان 
در جنوب‌شرق کشور بیشترین درصد افزايش مقدار بارش فرین پیش‌نگری شده نسبت به دوره تاریخی 
مشاهده می‌شود. در مطالعه زرین و داداشی رودباری ey (Zarrin & Dadashi-Roudbari, 2022a)‏ افزایش 
روزهای همراه با بارش سنگین در سواحل دریای خزر عنوان شده است و همچنین زرین و داداشی رودباری 
(Zarrin & Dadashi-Roudbari, 2022b)‏ نشان دادند که بيشینه بارش محتمل در دوره‌های بازگشت ۱ تا ۱۵ 
pelea asia‏ شمالی. lel pts‏ می ol,‏ بروسی dS‏ ۶ ایسگاه Ss‏ فده dabo du Ca. (St)‏ که 
بیشترین شیب درصد افزایش مقدار بارش فرین پیش‌نگری شده نسبت به مقدار بارش فرین در دوره تاریخی 
مربوط به ایستگاه‌های زاهدان و اراک است که در دوره باز CAS‏ ۲۰۰ ساله این درصد افزایش به بیش از ۷۰ 
درصد می‌رسد. در سناریو میانه (شکل CF‏ نسبت به دو سناریو دیگر در حالت کلی درصد افزایش مقدار 
بارش فرین پیش‌نگری شده نسبت به مقدار بارش فرین در دوره تاریخی بیشتر است. بیشترین مقادیر افزایش 
در رشته‌کوه زاگرس و سواحل جنوبی کشور (سنندج. همدان اراک» LEW‏ اصفهان» شیراز فساء میناب؛ 
بوشهر بندرعباس, بندرلنگه ابوموسی) و گرگان (شمال‌شرق کشور) واقع هستند. بر اساس شکل A)‏ 
بیشترین شیب درصد افزایش مقدار بارش فرین پیش‌نگری شده نسبت به مقدار بارش فرین در دوره تاریخی 
مربوط به ایستگاه‌های ابوموسی و اراک است که در دوره بازگشت ۲۰۰ ساله این درصد افزایش به‌ترتیب به 


حدود ۸۰ و ۶۸ درصد می‌رسد. در ایستگاه‌های جنوبی کشور (بندرعباس, بندرلنگه» بوشهر و میناب) از دوره 
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بازگشت ۱۰ سال به بعد شیب درصد افزایش مقدار بارش فرین پیش‌نگری شده در دوره بازگشت‌های 
طولانی‌تر کاهش می‌یابد. همانطور که در شکل Y‏ مشخص است کمترین در صد افزايش مقدار بارش فرین 
پیش‌نگری شده نسبت به مقدار بارش فرین در دوره تاریخی» در سناریو بدبینانه (SSP5-8.5)‏ حاصل شده 
است. بیشینه این افزايش مقدار بارش فرین در دوره آینده (۲۰۵۰-۲۰۲۶) مربوط به ایستگاه‌های بجنورد؛ 
بیرجند و ایستگاه‌های اراک DUIS‏ و سمنان واقع در مناطق مرکزی کشور است. در این سناریو شیب 
تغییرات JRE)‏ ۵6) نسبت به سایر سناریوها نیز کمتر است و بیشترین افزایش مربوط به ایستگاه اراک است 


که در دوره باز EE dus‏ ساله درصد افزایش مقدار بارش فرین نسبت به دوره تاریخی به حدود ۰ درصد 


می رسد. 
2year return period 55۳245: 2year return period SSP585: 2year return period‏ :55۴126 
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55۳126: 200۷90۲ return period 


55۳245: 200۷66۲ return period 


45E 505 55E 606 * 


45E 50 55E 60 * 


SSP585: 200۷90۲ return period 


شکل ۴- درصد تغییرات مقدار بارش فرین پیش‌نگری شده نسبت به مقدار بارش فرین در دوره تاریخی در هر 
یک از ایستگاه‌های موردمطالعه در سناریوهای SSP2-4.5 SSP1-2.6‏ 3 9۳5-5 و دوره با زگشت‌های at‏ 


۵ ۱۰ ۰۲۵ ۵۰ ۱۰۰ و ۲۰۰ Sle‏ ایستگاه‌های مشخص شده با خط بسته قرمزرنگ دارای درصد افزایش بیش از 


۲۵ درصد هستنك. 


Fig. 4. Percentage changes in the projected extreme precipitation amounts compared to 


historical extreme precipitation at each of the studied stations under the SSP1-2.6, SSP2-4.5, 
and SSP5-8.5 scenarios for return periods of 2, 5, 10, 25, 50, 100, and 200 years. Stations 
marked with a solid red line indicate a percentage increase of more than 25%. 
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شکل ۵- نمودار تغییرات کمیت ترسیم شده در شکل Y‏ در دوره بازگشت‌های ۲ ۸۵ ۰۱۰ ۲۵ ۵۰ ۱۰۰ و ۲۰۰ 
سال و در سنار یو SSP5-8.5 (c SSP2-4.5 (b SSP1-2.6 (a‏ برای ایستگاه‌های مشخص شده در شکل Y‏ 


Fig. 5. Chart of changes in the quantity plotted in Figure 4 for return periods of 2, 5, 10, 25, 50, 


100, and 200 years under scenarios a) SSP1-2.6, b) SSP2-4.5, and c) SSP5-8.5 for the stations 
indicated in Figure 4. 


۴-نتیجه گیری 


در پژوهش حاضر دوره بازگشت و احتمال رخداد بارش‌های فرین تا اواسط قرن E SL aS FLY‏ مدل 
مورداستفاده از مدل‌های CMIP6‏ و تحت سناریوهای SSP2-4.5 SSP1-2.6‏ و SSP5-8.5‏ موردبررسی و 
تحلیل قرار گرفت. با استفاده از آزمون نکویی برازش کای اسکوثر تابع توزیع احتمالاتی بیشینه بارش ۲۴ 
ساعته در دوره‌ی تاریخی (۲۰۱۴-۱۹۹۰) و دوره آینده (۲۰۵۰-۲۰۲۶) تحت سناریوهای تغییر اقلیم خوش - 


بینانه» میانه و بدبینانه (بر اساس کمترین خطای آزمون» بهترین توزیع Goll‏ منطبق بر داده‌ها برای هر یک از 


۳۳۲ نشریه جغرافیا و مخاطرات محیطی, جلد ۰۱۳ شماره ۲ پائیز ۱۳۰۳ 


ایستگاه‌های موردمطالعه انتخاب شد که در AF‏ درصد ایستگاه‌ها توزیع 3 log-Pearson‏ به‌عنوان بهترین 
ab‏ توزیع انتخاب گردید. نتایج بررسی بارش‌های فرین حاصل در هر ایستگاه و مقایسه آن‌ها در دوره آینده 
و تاریخی OLE‏ می‌دهد درصد تغییرات مقدار بارش فرین پیش‌نگری شده نسبت به مقدار بارش فرین در 
دوره تاریخی. در سناریو SSP2-45‏ نسبت به دو سناریو خوش‌بینانه و بدبینانه بیشتر است و در سناریو 
بدبینانه نست به ple‏ سناریوها کمترین درصد تغییرات مشاهده می‌شود. همچنین در سناریو خوش‌بینانه 
(SSPI-2.6)‏ بیشترین درصد افزايش مقدار بارش فرین پیش‌نگری شده نسبت به دوره تاریخی در ایستگاه 
گرگان و بابلسر واقع در سواحل شرقی دریای خزر و ایستگاه بجنورد در شمال‌شرق کشور و ایستگاه‌های 
اراک و اصفهان واقع در مرکز کشور و همچنین ایستگاه زاهدان در جنوب‌شرق کشور مشاهده می‌شود که 
نتایج بیشینه مقدار بارش فرین در سواحل دریای خزر در تطابق با مطالعات زرین و داداشی رودباری ( Zarrin‏ 
Dadashi-Roudbari, 20222: 0‏ &( است. در مطالعه زرین و داداشی رودباری Zarrin & Dadashi-)‏ 
acy (Roudbari, 2022b‏ شدت بارش در دوره‌های باز گشت ۱ OF‏ ۱۰ و ۱۵ ساله در سواحل جنوبی 
ois‏ زاگرس مرتفع و سواحل دریای خزر عنوان شده است؛ که در مقایسه با مطالعه حاضر بیشینه افزايش 
بارش فرین در سواحل جنوبی کشور در دوره بازگشت‌های طولانی‌تر مشاهده می‌شود. در سناریو خوش‌بینانه 
(SSPI1-2.6)‏ بیشترین شیب درصد افزایش مقدار بارش فرین پیش‌نگری شده نسبت به مقدار بارش فرین در 
دوره تاریخی مربوط به ایستگاه‌های زاهدان و اراک است که در دوره ناز کشت ۰ ساله این درصد افزایش 
به بیش از ۷۰ درصد می‌رسد. در سناریو SSP2-4.5) ale‏ بیشترین مقادیر افزایش در ایستگاه‌های سنندج» 
همدان. اراک کاشان. اصفهان. شیراز» فسا و برخحی مناطق سواحل جنوبی کشور شامل ایستگاه‌های میناب؛ 
Gedy‏ بندرعباس, بندرلنگه. ابوموسی و نیز در گرگان واقع در شمال‌شرق کشور مشاهده می‌شود. همچنین 
در این در سناریو, بیشترین شیب درصد افزايش مقدار بارش فرین پیش‌نگری شده نسبت به مقدار بارش 
فرین در دوره تاریخی مربوط به ایستگاه‌های ابوموسی و اراک است که در دوره بازگشت ۲۰۰ ساله این 
درصد افزایش به ترتیب به حدود ۸۰ و ۶۸ درصد می‌رسد. در سناریو میانه (SSP2-4.5)‏ در ایستگاه‌های 
سواحل جنوبی کشور شامل بندرعباس. بندرلنگه بوشهر و میناب از دوره بازگشت ۱۰ سال به بعد شیب 


درصد افزایش مقدار بارش فرین rye ued ego‏ در دوره b‏ گشت‌های طولانی تر کاهش LL‏ در سناریو 
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بدبینانه (85۳5-8.5) بیشینه افزایش مقدار بارش فرین مربوط به ایستگاه‌های بجنورد» بیرجند» اراک» کاشان و 
Obss‏ است. در این pile‏ شیب تغییرات نسبت به سایر سناریوها نیز کمتر است و بیشترین افزایش bpp‏ 
به ایستگاه اراک است که در دوره بازگشت ۲۰۰ ساله درصد افزایش مقدار بارش فرین نسبت به دوره تاریخی 
به حدود ۵۰ درصد می‌رسد. همچنین نتایج این مطالعه نشان می‌دهد در بین ایستگاه‌های موردبررسی افزایش 
مقدار بارش فرین نسبت به دوره تاریخی با درصد افزايش بیش از YO‏ درصد در ایستگاه اراک در هر سه 
سناریو. در ایستگاه‌های اصفهان و گرگان در سناریوهای SSP1-2.6‏ و ۹5۳2-45 و در ایستگاه بجنورد در 


سناریوهای SSP1-2.6‏ و SSP5-8.5‏ و کاشان در سناریوهای SSP2-4.5‏ و SSP5-8.5‏ تأیید شده است. 


در پژوهش حاضر دوره باز گشت و احتمال رخداد بارش‌های فرین در دوره ۲۰۵۰-۲۰۲۶ با ترکیب وزنی ۱۰ 
مدل مورداستفاده از مدل‌های 0111۳6 با در نظر گرفتن مهارت‌های شبیه‌سازی مدل در هر دو الگوی فضایی و 
تنوع بین‌سالانه و تحت سناریوهای SSP1-2.6‏ 85۳2-45 و 55۳5-85 موردبررسی و تحلیل قرار گرفت. با 
توجه به خطرات رویدادهای بارش فرین ممکن است در آینده این بارش‌های فرین در اکثر مناطق موردمطالعه 
به‌طور قابل‌توجهی افزایش یابد؛ بنابراین لازم است تلاش‌های بیشتری در GEIS‏ آینده برای بهبود روش 
وزن‌دهی و انتخاب مدل‌های بیشتر انجام شود. همچنین با توجه به دستاورد این پژوهش لازم است تا در هر 
منطقه تعداد بیشتری ایستگاه موردبررسی قرار گیرد بهخصوص در مناطقی مانند استان مرکزی که در این 
منطقه هر سه سناریو افزايش مقدار بارش فرین را نسبت به دوره تاریخی در دوره بازگشت‌های مختلف نشان 


داد. 
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